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ABSTRACT 

We present deep, high-quality iiT-band images of complete subsamples of powerful radio and 
sub-millimetre galaxies at redshifts z ~ 2. The data were obtained in the very best available 
seeing via queue-based observations at the United Kingdom InfraRed Telescope and Gemini 
North, with integration times scaled to ensure that comparable rest-frame surface brightness 
levels are reached for all of the galaxies. We fit two-dimensional axi-symmetric galaxy mod- 
els to these images to determine basic galaxy morphological parameters at rest-frame optical 
wavelengths Xrest > 4000A, varying luminosity, axial ratio, half-light radius ri/2, and Sersic 
index n. We find that, while a minority of the images show some evidence of galaxy inter- 
actions, > 95% of the rest-frame optical light in all the galaxies is well-described by these 
simple axi-symmetric models. We also find evidence for a clear difference in morphology be- 
tween these two classes of galaxy; fits to the individual images and to the image stacks reveal 
that the radio galaxies are moderately large {{ri/2) — 8.4 ± 1.1 kpc; median ri/2 = 7.8 ), 
de Vaucouleurs spheroids ((n) = 4.07 ± 0.27; median n = 3.87), while the sub-millimetre 
galaxies appear to be moderately compact {{ri/2) — 3.4 ± O.Skpc; median ri/2 = 3.1 kpc) 
exponential discs ((n) = 1.44 ± 0.16; median n = 1.08). Model fits to the stacked images 
yield very similar results. We show that the z 2± 2 radio galaxies display a well-defined Ko- 
rmendy (fig — Ve) relation but that, while larger than other recently-studied massive galaxy 
populations at comparable redshifts, they are still a factor ~ 1.5 times smaller than their local 
counterparts. The scalelengths of the starlight in the sub-millimetre galaxies are comparable 
to those reported for the molecular gas, suggesting that the two may be co-located. Their sizes 
are also similar to those of comparably massive quiescent galaxies at z > 1.5, allowing the 
possibility of an evolutionary connection following cessation/quenching of the observed star- 
formation activity. In terms of stellar mass surface density, the majority of the radio galaxies 
lie within the locus defined by local ellipticals of comparable stellar mass. In contrast, while 
best modelled as discs at the epoch of bright dust-enshrouded star formation, most of the sub- 
millimetre galaxies have higher stellar mass densities than local galaxies, and appear destined 
to evolve into present-day massive ellipticals. 

Key words: galaxies: active - galaxies: starburst - galaxies: photometry - galaxies: funda- 
mental parameters - infrared: galaxies. 



1 INTRODUCTION Driven by multi-wavelength surveys of ever increasing scale, 

recent years have seen substantial advances in our understanding 



The statistical and individual properties of galaxies selected at dif- 
ferent wavelengths, and at different redshifts must be explained in 
a self-consistent manner by any successful model of galaxy forma- 
tion and evolution. 



of the basic demographics of different galaxy populations such as 
Lyman-break galaxies (e.g. Bouwens et al. 2007), radio galaxies 
(e.g. Rigby et al. 2010) and sub-millimetre galaxies (e.g. Coppin et 
al. 2006; Weiss et al. 2009; Austermann et al. 2010). 



The resulting improving information in, for example, number 
* Email: tat@roe.ac.uk counts and redshift distributions, has been fed back to modify/refine 

t Scottish Universities Physics Alliance models of galaxy formation. However, important uncertainties and 
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degeneracies in tiie basic physics generally remain. For example, 
sub-millimetre galaxies in particular were not a natural prediction 
of early semi-analytic models (Baugh et al. 2005), and there re- 
mains controversy over whether this can be fixed by the adoption 
of a top-heavy stellar intial mass function (IMF), or whether this 
tension is symptomatic of a basic inability of current semi-analytic 
models to produce enough star-formation activity in the young uni- 
verse (Bower et al. 2006; Swinbank et al. 2008; Cirasuolo & Dun- 
lop 2008). 

To remove such uncertainties requires higher-level informa- 
tion. Specifically, different galaxy formation models make rather 
different predictions for galaxy masses, sizes and morphologies. 
For example, the models of Baugh et al. (2005) predict that sub- 
millimetre galaxies result from mergers of sub-L* galaxies (with 
their apparent star- formation rates boosted by a top-heavy IMF). 
In contrast, the recent cosmological hydrodynamic simulations of 
Dave et al. (2010) predict that sub-millimetre galaxies are massive 
galaxies fed by smooth infall, are not typically involved in major 
mergers, and are simply the most extreme examples of normal star- 
forming galaxies (see also Finlator et al. 2006; Fardal et al. 2007; 
Dekel et al. 2009). 

The estimation of galaxy stellar masses at high reshift has 
therefore become a major industry in recent years, boosted by 
the provision of rest-frame near-infrared photometry provided by 
the Spitzer Space Telescope. The sensitivity of mass-to-light ratio 
on star-formation history is minimized at such wavelengths, but 
uncertainties due to choice of IMF or the relative importance of 
asymptotic giant branch stars have motivated attempts to obtain in- 
dependent confirmation of galaxy masses at 2 ~ 2 via velocity 
dispersion measurements derived from emission-line spectroscopy 
at both millimetre and infrared wavelengths (e.g. Taconni et al. 
2006, 2008; van Dokkum, Kriek & Franx 2009). While acknowl- 
edging the uncertainties, the available photometric and dynamical 
evidence suggests that, by the epoch of peak sub-millimetre emis- 
sion at 2; ~ 2 — 3, sub-millimetre galaxies were already massive 
galaxies with estimated stellar masses Mstar — 10^^ M0 (Greve 
et al. 2004; Tacconi et al. 2006; Borys et al. 2005; Clements et al 
2008; Dye et al. 2008; Schael et al. 2010). This suggests that sub- 
millimetre galaxies represent an important phase in the formation 
of massive galaxies, and the challenge now is to establish their role 
in cosmic star-formation history, and the nature of their present- 
day descendents. Here we attempt to move this subject forward 
by measuring the basic sizes and morphologies of sub-millimetre 
galaxies at 2 ~ 2, performing such measurements at the longest 
possible wavelength (2.2 pm) in order to trace the mass dominant 
stellar population, and to minimize the potentially confusing effects 
of dust obscuration at rest-frame ultraviolet wavelengths. 

The measurement of galaxy sizes and morphologies at red- 
shifts 2 > 1 has also been the subject of many studies in recent 
years, and has yielded some surprising and still controversial re- 
sults. Specifically, a number of studies have now found that mas- 
sive, passive galaxies at 1.5 < 2 < 3 are much more compact 
than present-day galaxies of comparable mass (Trujillo et al. 2006, 
2007; Longhetti et al. 2007; Cimatti et al. 2008; van Dokkum et al. 
2008). The resulting inferred high stellar mass densities have been 
taken to imply that such objects cannot continue to evolve by purely 
passive evolution into any known type of present-day galaxy. These 
results are consistent with those of Toft et al. (2007) and Zirm et 
al. (2007) who, with the aid of Spitzer MIPS photometry, were 
in addition able to show that, while passive high-redshift galax- 
ies are surprisingly compact, star-forming galaxies at high redshift 
appear morphologically similar to normal star-forming disk galax- 



ies. Recently Buitrago et al. (2008), from a Hubble Space Tele- 
scope (HST) NICMOS study of 82 massive galaxies in GOODS 
North concluded that, for a given stellar mass, disc-like galaxies 
at 2 ~ 2.3 were a factor 2.6 ± 0.3 smaller than present-day sys- 
tems, while spheroid-like galaxies were a factor 4.3 ± 0.7 smaller. 
They also noted that this evolution is at least in qualitative agree- 
ment with some hierarchical semi-analytic predictions (e.g. Khoch- 
far & Silk 2006). However, there is now growing evidence that not 
all high-mass galaxies at 2 > 1 follow this trend, with some at 
least appearing to be fully formed (Saracco, Longhetti & Gargiulo 
2010). 

At present it is unclear where sub-millimetre galaxies fit into 
this picture. As star-forming galaxies they might be expected to 
be extended discs, like Lyman-break galaxies and present-day star- 
forming galaxies. On the other hand, if they are the progenitors 
of today's massive ellipticals, their apparently high-density star- 
bursts (Tacconi et al. 2006) might be related to the high stellar den- 
sities reported in red, passive galaxies at high redshift (~ 2 x 
lO^^Mo kpc^^, comparable to present-day globular clusters). 

To date, genuinely deep, if-band imaging has only been ob- 
tained for a small number of sub-millimetre galaxies (e.g. Small 
et al. 2003), and generally with inadequate resolution for the reli- 
able extraction of basic morphological parameters. Very recently, 
Swinbank et al. (2010) reported the results of HST NICMOS-^ACS 
imaging of a sub-sample of sub-millimetre galaxies. However, they 
did not report reliable measurement of, for example Sesic indices 
for individual galaxies, and the wavelength limit of HST prevented 
them from undertaking observations in the ii'-band, which ensures 
imaging longward of the 4000A/Balmer break out to 2 ~ 3. 

A high-quality, systematic, deep, if-band imaging study of a 
well defined sample of sub-millimetre galaxies is thus desirable, 
and has the potential to clarify several fundamental issues. In this 
paper we attempt to take a step towards better quantifying the sizes, 
morphologies and stellar mass densities of sub-millimetre galax- 
ies via a deep, and high-quality Tf-band imaging study of a well- 
defined sample of luminous SMGs at 2 ~ 2. To aid the interpre- 
tation of our results we have taken the somewhat unusual step of 
obtaining images of equivalent depth and resolution for a complete 
comparison sample of the most luminous radio galaxies known at 
comparable redshifts (2 ~ 2). 

The radio galaxy sample is a useful control sample because, 
given the now-well-established proportionality of black hole and 
spheroid mass (e.g. McLure & Dunlop 2002), we can be very con- 
fident that these galaxies will end up as massive ellipticals at the 
present-day. However, given the above-mentioned evidence that 
many massive galaxies at 2 > 1 are more compact than compa- 
rably massive present-day objects, a study of radio galaxy sizes at 
2 ~ 2 is also clearly of interest in its own right. Rather surprisingly, 
while deep near-infrared imaging of 3CR radio galaxies has been 
undertaken at 2 ~ 1 (e.g. Dunlop & Peacock 1993; Best et al. 1998) 
very deep, high-quality ii'-band imaging of the most powerful 3CR 
radio galaxies at 2 ~ 2 has not in fact been undertaken before (al- 
though NICMOS i?-band images of a sample of even higher red- 
shift radio galaxies were obtained by Penterrici et al. 2001). Conse- 
quently, while the morphologies of radio galaxies out to 2 ~ 1 are 
rather well defined (McLure & Dunlop 2000; McLure et al. 2004), 
their rest-frame optical morphologies and sizes at 2 ~ 2 have yet 
to be clarified. 

To undertake this study we have used NIRI on Gemini North 
to image the sub-millimetre galaxies, and UFTI on UKIRT to im- 
age the radio galaxies. We have taken great care to ensure that 
all images are as well matched as possible in K-hanA surface- 
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brightness limit, designing integration times to compensate for ttie 
difference in telescope aperture between Gemini North and UKIRT, 
and for cosmological surface-brightness dimming within our sam- 
ples. Also crucial to this study has been our ability (due to the 
queue-based scheduling system in operation at both telescopes) to 
demand the very best seeing conditions for all of our observations 
(FWHM < 0.6 arcsec). It is only with this powerful combination 
of consistent high-quality seeing and long exposure times (typi- 
cally 90 mins on Gemini) that it has proved possible to determine 
the sizes (i.e. half-light radii) and basic morphologies (i.e. Sersic 
index) of these massive z ~ 2 galaxies via ground-based imaging. 
The paper is structured as follows. In Section 2 we summarize 
the radio galaxy and sub-millimetre galaxy samples, and the exist- 
ing supporting information. In Section 3 we present the UKIRT and 
Gemini if-band imaging observations, and in Section 4 we explain 
how these data were reduced. Then in Section 5 we describe how 
we modelled these new images to extract the basic morphological 
properties of the galaxies. The results of the individual model fits 
are presented, summarized and discussed in Section 6, with the im- 
ages, model fits and residual images provided in Appendix A. We 
also explore the robustness of our results, through simulations, non 
parametric tests of galaxy size, and the analysis of radio galaxy 
and sub-millimetre galaxy image stacks. Finally, in Section 7 we 
attempt to place our results in the context of studies of other galaxy 
populations at both high and low redshift. Our main conculsions 
are summarized in Section 8. 



2 SAMPLES AND SUPPORTING DATA 

The initial target sample for this study consisted of 39 sources, and 
comprised a complete sample of the 14 most radio-luminous ra- 
dio galaxies at 2: ~ 2, and the 25 brightest 850 /im sources in the 
SCUBA 8-mIy survey. 



2.1 The radio galaxy sample 

Of the 14 objects in the original radio galaxy sample, 8 were drawn 
from the 3CRR catalogue (Laing, Riley & Longair 1983) and 6 
from the equatorial sample of Best et al. (1999). The 3CRR sample 
is complete to a 178-MHz flux-density limit of Sits mhz > 10.9 Jy 
over an area of 4.2 steradians, and contains a total of 170 steep- 
spectrum sources over the redshift range 0.0 < z < 2.5. These data 
are essentially doubled by the addition of the Best et al. sample, 
which was designed to have a roughly equivalent selection crite- 
rion of 5*408 MHz > 5 ly, using the Molonglo Reference Catalogue. 
From this combined survey we have defined a complete (cf Pen- 
tericci et al. 2001) sample of the 14 most powerful, low-frequency- 
selected, radio sources in the redshift range 1.5 < z < 2.5. One of 
the sources in this high-power sample, 3C 294, transpired to lie too 
close on the sky to a nearby star for deep if-band imaging to be 
successfully secured. Thus, for the present study, the parent sample 
was trimmed to the 13 sources listed in Table 1. 



2.2 The sub-millimetre galaxy sample 

The sub-millimetre galaxies were drawn from the SCUBA 8-mIy 
survey (Scott et al. 2002). The SCUBA 850 ^m imaging in the 8- 
mjy survey was split between the Lockman Hole East and Elais 
N2 fields, and covered a total area of 260 arcmin^ to a typical 
rms noise level of asso — 2.5mJybeam~^. From the resulting 



source catalogue, all 25 sub-millimetre sources with a signal-to- 
noise S/N > 3.5 at 850 ^im were selected as potential targets for 
deep if-band imaging in the present study. In Table 1 we list the 15 
sub-millimetre sources which were actually observed, and detected 
with sufficient signahnoise in if-band to undertake a meaningful 
morphological analysis. 

2.3 Identification of K-band counterparts 

All of the radio galaxies in the 13-source sample have unambigu- 
ous near-infrared galaxy counterparts, as a result of the sub-arcsec 
positional accuracy provided by the VLA detections of their radio 
cores. 

However, as is well known, the identification of the A'-band 
counterparts of the sub-millimetre sources is more complicated, 
due to the large (FWHM ~ 15 arcsec) size of the Airy disc de- 
livered by the 15-m JCMT at A ~ 850 /im. However, for 12 of the 
15 sub-millimetre sources listed in Table 1, robust radio identifica- 
tions (and hence positions accurate to ~ 1 arcsec) were successfully 
deduced by Ivison et al. (2002), via deep 1.4 GHz imaging with the 
Very Large Array (VLA) reaching a typical 5a detection limit of 
SiAGHz — 25/itJy. 

As indicated in Table 1, identifications for the remaining 3 
sub-millimetre sources in our sub-sample were therefore by neces- 
sity determined on the basis of optical-infrared information only. 
For the 2 outstanding ELAIS N2 sources, the choice of identifica- 
tion was based on the I — K colour derived from the existing (and 
now public) /-band HST ACS imaging of the ELAIS N2 field (PI: 
Almaini, PID: 9761), and the pre-existing UKIRT A'-band imaging 
described by Ivison et al. (2002). Specifically, following Pope et al. 
(2005) the criterion that sub-millimetre galaxies are typically very 
red in I — K colour (see also Small et al. 2004) was used to se- 
lect the most likely K-hanA counterpart in the absence of a secure 
radio detection. In the case of LE 850.04, where no optical imag- 
ing was available at the commencement of this study, the ii'-band 
source closest to the sub-millimetre position was selected. It should 
be noted that although LE 850.04 is therefore potentially the most 
uninformed identification, no other possible A'-band counterpart is 
found within the SCUBA positional error circle for this source. 

While it is to some extent desirable to partially remove the 
biases introduced by insisting on a robust radio identification, it is 
nonetheless true that these three non-radio selected galaxy identi- 
fications could be regarded as less secure than their radio-selected 
counterparts. Because of this, they are explicitly flagged as a spe- 
cific subset in much of the subsequent analysis presented in this 
paper. 

2.4 Redshifts 

Secure spectroscopic redshifts exist for all 13 radio galaxies (gen- 
erally aided by bright emission lines) and are given in Table 1. 

As is well documented in the literature (e.g. Ivison et al. 2005), 
obtaining secure spectroscopic redshifts for sub-millimetre galax- 
ies has proved much more difficult, not only because the identifi- 
cations are harder to establish, but also because the correct galaxy 
identifications are often very red, providing little optical light for 
optical spectroscopy. Where secure identifications exist, spectro- 
scopic redshifts are thus typically only available for sub-millimetre 
sources with optically bright counterparts or strong-lined AGN 
emission (Chapman et al. 2005). Nevertheless secure spectroscopic 
redshifts have now been secured for 8 of the 15 sub-millimetre 
sources, and these are listed in Table 1. 
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Table 1. Details of the radio and sub-millimetre galaxy _R'-band observations. Column 1 lists the lAU name for the radio galaxies, and the 8-mJy survey 
name for the sub-millimetre sources. Column 2 gives the original Cambridge and Parkes catalogue names for the radio sources, and the alternative SHADES 
catalogue names for the Lockman sub-millimetre galaxies (Coppin et al. 2006). Column 3 provides information on the nature of the galaxy identification 
for the sub-millimetre sources. Column 4 names the telescope with which the data were acquired, and Column 5 lists the total exposure time in seconds for 
each target. Column 6 provides the measured 6-arcsec diameter aperture K-hanA magnitude (in the Vega magnitude system) for each galaxy, along with the 
associated error Finally, Columns 7, 8, and 9 list the adopted redshift for each source, the type of redshift, and the reference for the source from which this 
redshift was obtained: A from Best el al. (1999), B from Hewitt & Burbidge (1991), C from Spinrad et al. (1985), a from Aretxaga et al. (2005), b from 
Aretxaga, priv. comm., c from Chapman et al. (2003), d from Chapman et al. (2005), e from Ivison et al. (2005), f from Small et al. (2003). 



Source 


Radio Catalogue Name 




Telescope 


Exp. time (sec) 


K mag. (Vega) 


z 


2-type 


z-ref 


0016-12 


3C 008 




UKIRT 


9480 


18.13+0.07 


1.589 


Zspec 


A 


0128--26 


PKS 0128-26 




UKIRT 


18840 


17.77+0.09 


2.348 


Zspec 


A 


0231+31 


3C 068.2 




UKIRT 


9720 


17.97+0.08 


1.575 


^spec 


B 


0310-15 


PKS 0310-15 




UKIRT 


8100 


18.35+0.11 


1.769 


Zspec 


A 


0851-14 


PKS 0851-14 




UKIRT 


10140 


17.87+0.08 


1.665 


^spec 


A 


1008+46 


3C239 




UKIRT 


11820 


17.89+0.07 


1.790 


Zspec 


C 


1019+22 


3C241 




UKIRT 


9600 


17.75+0.08 


1.617 


Zspec 


C 


1120+05 


3C257 




UKIRT 


22560 


17.58+0.07 


2.474 


Zspec 


B 


1140-11 


PKS 1140-11 




UKIRT 


13920 


18.79+0.09 


1.935 


Zspec 


A 


1422-29 


PKS 1422-29 




UKIRT 


9600 


17.69+0.07 


1.632 


Z'spec 


A 


1533+55 


3C322 




UKIRT 


9180 


17.86+0.08 


1.681 


Z'spec 


B 


1602-17 


PKS 1602-17 




UKIRT 


15120 


17.63+0.09 


2.043 


^spec 


A 


2356+43 


3C 470 




UKIRT 


10260 


18.11+0.12 


1.653 


^spec 


B 


Source 


SHADES Name 


ID Type 


Telescope 


Exp. time (sec) 


K mag. (Vega) 


Z 


z-type 


z-ref 


N2 850.01 




Radio ID 


Gemini N 


4750 


19.90+0.15 


2.850 


Zspec 


d 


N2 850.02 




Radio ID 


Gemini N 


4500 


20.77+0.27 


2.454 


Zspec 


d 


N2 850.04 




Radio ID 


Gemini N 


4625 


18.85+0.10 


2.387 


Zspec 


f 


N2 850.06 




EROA'RO 


Gemini N 


5500 


20.18+0.16 


3.0 


^phot 


b 


N2 850.07 




Radio ID 


Gemini N 


5000 


20.00+0.16 


3.1 


Zphot 


b 


N2 850.08 




Radio ID 


Gemini N 


5000 


18.30+0.09 


1.190 


Zspec 


d 


N2 850.09 




Radio ID 


Gemini N 


4750 


20.99+0.37 


3.3 


^phot 


b 


N2 850.12 




EROA'RO 


Gemini N 


5500 


20.31+0.17 


2.5 


^phot 


a 


N2 850.13 




Radio ID 


Gemini N 


5375 


21.61+0.37 


2.283 


Z'spec 


d 


LE 850.02 


LOCK850.3 


Radio ID 


Gemini N 


5500 


19.66+0.11 


3.041 


Z'spec 


d 


LE 850.03 


LOCK850.17 


Radio ID 


Gemini N 


5500 


19.66+0.12 


2.239 


Z'spec 


d/e 


LE 850.04 


LOCK850.6 


EROA'RO 


Gemini N 


5500 


20.74+0.14 


3.6 


^phot 


b 


LE 850.06 


LOCK850.14 


Radio ID 


Gemini N 


5625 


19.16+0.17 


2.610 


Zspec 


c/e 


LE 850.07 


LOCK850.16 


Radio ID 


Gemini N 


5500 


18.87+0.10 


1.9 


Zphot 


b 


LE 850.08 


LOCK850.41 


Radio ID 


Gemini N 


5500 


20.38+0.16 


3.4 


^phot 


b 



In the absence of secure spectroscopic redshifts, photometric 
redshifts based on the observed shape of the sub-millimetre-radio 
spectral energy distribution have been adopted from Aretxaga et al. 
(2007). These redshift estimates have been found to be accurate to 
Az < 0.5 when compared to existing spectroscopic values. We 
note that in the general high-redshift range of interest here, derived 
physical galaxy scalelength is relatively insensitive to the precise 
redshift. 



that even passive evolution of the stellar population offsets to some 
extent the impact of surface brightness dimming with increasing 
redshift. In the design of the observations, galaxy scalelengths of 
lOkpc were assumed (comparable to those found by McLure & 
Dunlop (2000) for 3CR galaxies atz = 1), and the total anticipated 
it -band magnitudes were estimated using the infrared Hubble dia- 
gram of De Breuck et al. (2002). 



3 OBSERVATIONS 

3.1 UKIRT observations of the radio galaxies 

The z ~ 2 radio galaxies were imaged in the K-band with the 3.8- 
m United Kingdom Infrared Telescope (UKIRT) on Mauna Kea. 
All observations used the UKIRT Fast-Track Imager (UFTI), a 1- 
2.5 fim camera with a 1024^ HgCdTe array, a pixel scale of 0.091 
arcsec, and hence a 1.5' x 1.5' field of view. All data were taken in 
< 0. 6-arcsec seeing and photometric conditions. Integration times 
were set to obtain a signal-to-noise of 5 at twice the anticipated 
half-light radius in the azimuthally binned luminosity profile. Note 



3.2 Gemini observations of the sub-mllllmetre galaxies 

The sub-millimetre sources were imaged in the K—band with the 
8-m Gemini North Telescope on Mauna Kea. All observations used 
the Near-Infrared Imager (NIRI), a 1-5 /im camera with a 1024^ 
ALADDIN InSb array. Camera f/6 was used, delivering a pixel 
scale of 0.1171 arcsec, and hence a 2.0' x 2.0' field of view. All 
data were taken in < 0.6 arcsec seeing and photometric condi- 
tions. Integration times were near constant over all observations at 
~ 1.5 hours per source, a choice made to yield images of compa- 
rable depth to those obtained with the longest exposures used on 
UKIRT for the most distant (z > 2) radio galaxies in the compari- 
son sample (see Table 1). The resulting it-band imaging reaches a 
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3(7 limiting magnitude of K 
aperture. 



4 DATA REDUCTION 



22.5 (Vega) in a 4-arcsec diameter 



Data reduction was performed using standard IRAF packages. Dark 
frames of equal integration time to the science data were taken on 
all nights. Each dark-subtracted image was divided by a normalized 
flat-field, derived though the scaled sigma-clipped median combi- 
nation of neighbouring science frames. Registering of the offset 
frames was performed using the brightest stars in the images. A 
map of the bad pixels in the array was obtained and used to exclude 
these pixels during image combination, followed by cosmic-ray re- 
jection. The initial reduced image was processed to create a mask 
of all source flux detected in the mosaiced data. With all known 
source flux excluded, the science frames were then reprocessed to 
create an improved flat-field, which was used in a re-reduction of 
the data to produce the final science frames. 



4.1 Astrometry and Photometry 

Given the numerous faint sources visible in these deep data, accu- 
rate astrometry was essential for associating the Ivison et al. (2002) 
radio positions with the correct iC'-band counterparts. Owing to un- 
certainties on scales of a few arcsec in the telescope pointing, it was 
therefore necessary to re-determine the astrometry of the images. 
Using GAIA, finding-charts from the Two Micron All Sky Survey 
(2MASS) were obtained, and the co-ordinates of stars present in 
both the 2MASS data and the K-h&nA imaging were used to deter- 
mine the plate solution for each image. After correction, the offsets 
between coordinates of 2MASS sources and their counterparts in 
the A'-band data were consistently < 0.1 arcsec. 

The aperture photometry listed in Table[T]was obtained using 
the IRAF APPHOT package. Background sky counts were mea- 
sured and subtracted using a curve-of-growth analysis. Faint stan- 
dard stars were observed for each target object for photometric cal- 
ibration. These data were reduced via the process applied to the sci- 
ence data. The observations of these standards consistently yielded 
zero-points of 22.3 ± 0.06 for UFTI and 23.68 ± 0.08 for NIRI, 
consistent with the values given in the literature. 



4.2 Point spread function 

A point spread function (PSF) for each science frame was built 
using the IRAF package PSF for use during two-dimensional mod- 
elling. Stars of comparable brightness to the target were selected 
and fitted with a range of alternative statistical distributions to de- 
termine the best fit. The resulting models were then centroided, av- 
eraged, and processed with look-up tables to reproduce any asym- 
metries in the PSF. 



5 HOST GALAXY ANALYSIS 

5.1 Image modelling 

A tried and tested two-dimensional modelling code (previously de- 
veloped for the more demanding task of studying quasar host galax- 
ies) was used to extract basic galaxy morphological parameters 
from the AT-band images of the sub-millimetre and radio galaxies. 
Full details of the two-dimensional modelling code can be found in 



McLure et al. (2000). In brief, the code fits the image with an ax- 
ially symmetric model (convolved with the observed PSF) with 5 
free parameters (luminosity, half-light radius, position angle, axial 
ratio, and Sersic index n), minimising x^ using a downhill-simplex 
technique. In order to ensure that the solution was stable (and not a 
local minimum), the modelling code was repeatedly restarted from 
a range of positions in the parameter space close to the minimum 

further below in Section 6. 



5.2 Reliability and limitations of modelling results 

The accuracy of scalelength and Sersic index recovery was exam- 
ined via two-dimensional modelling of synthetic galaxies. Model 
galaxies at 2 ~ 2 were constructed with a scalelength ri/2 of 1, 
1.5, 2, 3, 4, 5, 6, 7, and 8 kpc, and with Sersic index values of n = 1 
(exponential disc) and n = 4 (de Vaucouleurs spheroid). The axial 
ratio and position angle were set to be constant, with luminosity and 
redshift adopted from the mean values of the radio-galaxy sample. 
The synthetic sources were inserted into an image, with noise lev- 
els based upon the average rms-per-pixel values derived from the 
real images. Background counts were also added to each simulated 
model using the IRAF routine MKNOISE. 

To ensure a realistic level of uncertainty in the accuracy of the 
PSF used to model a given galaxy, synthetic data were convolved 
with a PSF derived from a science image, and then the modelling 
was undertaken using a new PSF constructed from a selection of 
stars elsewhere in the same field. 

Example results of these tests are given in Tables [3] and |4] 
These tests were repeated using PSFs and sky noise from other tar- 
gets in the UKIRT and Gemini samples with similar results, demon- 
strating that host galaxy scalelength and Sersic index can be recov- 
ered to within a typical accuracy of 10 per cent. It was, however, 
also found that Sersic index can only be usefully recovered when 
ri/2 > 2 kpc; for smaller galaxies the limitations of ground-based 
seeing render ellipticals and disc galaxies virtually indistinguish- 
able. 



6 RESULTS 

6.1 Scalengths and luminosity profiles 

The original iiT-band images, the resulting best fitting axi- 
symmetric model, and the residual images are presented in Ap- 
pendix A (Fig. Al for the radio galaxies, and Fig. A2 for the sub- 
millimetre galaxies). The best-fitting values for the model parame- 
ters are summarized in Table 2. 

A rather clear separation is found in the sizes and luminos- 
ity profiles of these two classes of galaxy, as can be seen from 
inspection of Table |2] The high-redshift radio galaxies emerge as 
moderately-large dynamically-relaxed ellipticals, while the sub- 
millimetre sources appear better described as disc galaxies. A com- 
parison of the Sersic indices is presented in Fig.[T] The 13/15 sub- 
millimetre galaxies with reliable modelling results (i.e. excluding 
N2 850.2 and N2 850.6 which have ri/2 < 2 kpc) possess Sersic 
values peaked at n = 1, and confined to the range 1 < n < 2.5. 
The Sersic indices of the 11/13 reliably-modelled radio galaxies 
(i.e. excluding 1 120+057 and 1602-174; see below) are in the range 
2.5 < n < 5.75 and peaked around n = 4, corresponding to a de 
Vaucouleurs luminosity profile. The distribution of half-light radii 
is shown in Fig.|2] Sub-millimetre galaxy scale-lengths are clearly 
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Table 2. Results from two-dimensional modelling of the ftT-band images of the radio galaxies (upper section) and the sub-millimetre galaxies (lower section). 
Column 1 gives the source name. Column 2 lists the semi-major axis scalelength (half-light radius) of the host galaxy fit in kiloparsec. Column 3 lists the mean 
surface brightness in K-hand magnitude (Vega) per square arcsec. Column 4 lists the value of the Sersic index in). Column 5 gives the axial ratio of the host 
galaxy. Column 6 lists the total host-integrated ii'-band magnitude (Vega). Column 7 gives the value of reduced x J for each model fit, and finally Column 8 
gives the % of K-hand flux from the object which is attributed to the smooth, axi-symmetric model. 



Source 


''1/2 


<Me > 


Sersic n 


Axial ratio 


K—nmg 


xl 


% flux in model 


0016-129 


7.58 


19.03 


5.23 


0.52 


17.84 


1.01 


98.3 


0128-264 


15.66 


20.64 


3.28 


0.45 


17.37 


0.94 


98.5 


0231+313 


11.54 


19.49 


2.73 


0.69 


17.70 


0.95 


99.3 


0310-150 


7.75 


19.51 


4.51 


0.33 


18.07 


0.82 


98.2 


0851-142 


9.72 


19.78 


3.85 


0.68 


18.19 


0.80 


96.9 


1008+467 


4.24 


17.44 


3.76 


0.85 


17.81 


0.94 


99.4 


1019+222 


5.09 


17.97 


4.59 


0.60 


17.72 


0.96 


97.7 


1120+057 


11.85 


16.69 


6.26 


0.42 


17.44 


1.14 


88.1 


1140-114 


4.20 


19.39 


3.05 


1.00 


18.90 


1.02 


95.6 


1422-297 


9.29 


19.16 


5.74 


0.45 


17.40 


0.97 


99.5 


1533+557 


5.95 


18.90 


4.21 


1.00 


17.82 


0.93 


98.9 


1602-174 


0.90 


14.99 


4.67 


0.88 


17.84 


0.88 


99.8 


2356+438 


11.02 


20.11 


3.87 


0.71 


18.47 


0.99 


97.3 


N2 850.1 


2.09 


18.61 


1.74 


0.76 


19.77 


0.99 


97.7 


N2 850.2 


1.15 


17.78 


3.33 


0.47 


20.83 


1.06 


95.9 


N2 850.4 


3.24 


18.2 


1.79 


0.66 


18.65 


0.85 


99.3 


N2 850.6 


1.95 


18.14 


3.29 


0.37 


20.19 


0.90 


96.1 


N2 850.7 


2.77 


19.16 


1.08 


0.75 


19.67 


0.72 


96.2 


N2 850.8 


3.70 


18.09 


1.35 


0.72 


18.21 


0.79 


99.4 


N2 850.9 


3.05 


19.74 


1.08 


0.41 


20.65 


0.79 


96.0 


N2 850.12 


4.65 


19.11 


2.17 


0.19 


20.06 


0.93 


99.5 


N2 850.13 


5.78 


20.34 


2.75 


0.22 


20.76 


0.79 


96.9 


LE 850.2 


3.88 


18.24 


1.07 


0.19 


19.48 


0.82 


97.2 


LE 850.3 


2.56 


17.53 


0.97 


0.53 


18.75 


0.58 


96.7 


LE 850.4 


2.87 


18.76 


0.93 


0.20 


20.49 


0.96 


96.1 


LE 850.6 


3.60 


19.3 


1.99 


0.89 


19.16 


0.90 


97.8 


LE 850.7 


3.10 


17.89 


0.86 


0.55 


18.71 


0.80 


97.7 


LE 850.8 


3.04 


20.04 


0.94 


0.75 


20.16 


0.72 


97.2 



peaked around r\/2 ~ 3kpc, while the radio galaxies are spread 
over a larger range, and are systematically larger (with an average 
'^1/2 = 8kpc). These results are discussed further in Section 7. 

6.2 Notes on individual sources 

6.2.1 Radio galaxies 

As explained in Section 2, all the z ~ 2 radio galaxy targets 
were initially identified with known radio positions in the litera- 
ture (Laing, Riley & Longair 1983, Best et al. 1999). In the light 
of our deep A'-band imaging, we have revised one of these radio- 
galaxy identifications. This is the case of 0851 — 14, where Best et 
al. (1999) selected a discy object to the east of the radio position 
as seen in their i?-band imaging. Our deep A"-band data reveals a 
multi-component merging system closer to the radio centroid. The 
larger component of this merger (which also displayed a redder 
R — K colour) was adopted as the true galaxy counterpart, and it 
is this object which is shown and modelled in Fig. Al. 

In the case of 1120+05, the highest-redshift radio galaxy in 
the 3CR catalogue, we find a partially-obscured point-like nucleus 
which is clearly shown in the two-dimensional modelling resid- 
ual (see Fig. lAlt . A nuclear component is also clearly present in 
1602—17. In this case, the i?-band image from Best et al. (1999) 
shows a distinctly different morphology (with no sign of activity in 
the nucleus) to that revealed in our i^-band data, which is domi- 



nated by a very compact source suggestive of obscured AGN emis- 
sion. As we do not believe our data are of sufficient quality to reli- 
ably separate a significant nuclear contribution from the host, and 
still derive robust galaxy scalelengths and Sersic indices, 1 120+05 
and 1602—17 have been excluded from further discussion of the 
statistical distribution of galaxy morphological parameters for the 
radio-galaxy sample. 



6.2.2 Sub-mm galaxies 

Due in part to their early redshift determinations by Chapman et 
al. (2003; 2005), two of the ELAIS N2 8-mJy sources, N2 850.2 
and N2 850.4, have previously been imaged in some detail at sub- 
millimetre, optical and near-infrared wavelengths (Small et al. 
2003; Tacconi et al. 2006; Tacconi et al. 2008). In particular, Tac- 
coni et al. (2008) report that, as viewed via high-resolution mm- 
wave CO line emission, N2 850.2 is a very compact object, with 
'"1/2 = 0.8 ± 0.5 kpc, while N2 850.4 is somewhat more extended, 
with riy2 = 2.4 ± 1 kpc. These figures agree remarkably well with 
the completely independent (and, at least in principle, physically 
distinct) measures of the half-light radius of the starlight in these 
two galaxies, as derived here from our modelling of the J-i'-band 
imaging (N2 850.2 has the smallest value of stellar half-light radius 
listed in Table 2, with ri/a = 1.1 kpc, while for N2 850.4 the de- 
rived value is ri/2 ~ 3.2 ± 0.3 kpc). This suggests that the basic 
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Table 3. Example results of two-dimensional modelling tests using syn- 
thetic disc galaxies (i.e. with a Sersic index n = 1). Column 1 lists the 
actual input scale-lengths of the simulated sources. Columns 2 and 3 list 
the recovered scale-lengths and Sersic indices as derived from the two- 
dimensional modelHng code, ri/2 values are given in parsecs at ^ = 2. 



Simulated rj ;2 


Recovered r^ ;2 


Recovered n 


1000 


878 


i.n 


1500 


1389 


1.97 


2000 


1886 


1.52 


3000 


3089 


1.12 


4000 


4181 


1.07 


5000 


5246 


1.04 


6000 


6375 


1.06 


7000 


7180 


1.02 


8000 


8238 


1.00 


9000 


9138 


1.10 


10000 


10192 


1.00 



Table 4. Example results of two-dimensional modelling tests using syn- 
thetic elliptical galaxies (i.e. with a Sersic index n = 4). Column 1 lists the 
actual scalelengths of the simulated sources. Columns 2 and 3 list the recov- 
ered scalelengths and Sersic indices as derived from the two-dimensional 
modelling code, ri/2 values are given in parsecs at 2 = 2. 



Simulated ri/2 


Recovered ri/2 


Recovered n 


1000 


882 


2.64 


1500 


1402 


3.29 


2000 


1920 


3.51 


3000 


3256 


3.89 


4000 


4115 


3.74 


5000 


5220 


3.92 


6000 


5936 


4.05 


7000 


6185 


3.94 


8000 


8223 


3.79 


9000 


9082 


3.82 


10000 


10076 


3.90 



distribution of the stars and molecular gas in these galaxies may be 
rather similar. 

Tacconi et al. (2008) also present a colour optical-infrared im- 
age for N2 850.4, derived from HST ACS and NICMOS imaging 
of this sub-mm galaxy. This shows a complex morphology, but as 
we demonstrate in sub-section 6.6, this is not inconsistent with our 
finding that the /("-band light is well-described by a rather straight- 
forward, axi-symmetric smooth galaxy model. Rather it provides 
further evidence of the dramatic wavelength dependence of mor- 
phology in at least some sub-mm galaxies (due to localized star- 
formation and/or patchy obscuration), and re-emphasizes the im- 
portance of moving to the longest possible near-infrared wave- 
length to uncover the basic morphological parameters of the mass- 
dominant host-galaxy stellar population. 



progressively larger apertures out to a radius of 12 arcsec. Upon de- 
termining the extent of the total flux from the source, the value of 
ri/2 was read from the profile. The Petrosian radius was defined as 
described by Blanton et al. (2001) and Yasuda et al. (2001), taken 
to be the point at which the ratio of the surface brightness in an an- 
nulus at radius r to the mean surface brightness within an aperture 
of radius r equals 0.2. 

The values of ri/2 as determined from the two-dimensional 
modelling, the curve-of-growth, and Petrosian techniques were 
found to be consistent with within 10 per cent. Typically, the results 
from two-dimensional modelling are systematically smaller than 
those derived from curve-of-growth techniques. This is expected 
given that the two-dimentional modeling effectively deconvolves 
the PSF (typically ~ 4kpc) from the data, yielding an unbiased 
measurement of true galaxy size. 



6.4 Non-parametric tests of galaxy morphology 

The distributions of the axial ratios and concentration indices for 
the two samples are plotted as a model-independent test of galaxy 
morphology. The axial ratios are plotted in Fig. [3] Although de- 
tailed statistical comparison is limited by small sample size, it can 
be seen that the distribution of axial ratios for the sub-millimetre 
galaxies is very broad and flat, extending to values ~ 0.1, as ex- 
pected for a population of disc galaxies (Sandage et al. 1970). In 
contrast, the axial ratios of all but one of the radio galaxies are con- 
fined to > 0.45, and the mean value of axial ratio is ~ 0.7, more 
consistent with that expected from a sample drawn from a popula- 
tion of elliptical galaxies (Sandage et al. 1970). 

The distribution of source concentration index is shown in Fig. 
[4] Concentration index represents the fraction of stars in the bulge 
component of a galaxy (Conselice 2003), and was calculated using 
Equation [T] where rgo and r2o are the radii containing 80 per cent 
and 20 per cent of the total light from the galaxy. 



C = 5 X log(r8o/r-2o) 



(1) 



The lower values found for the sub-millimetre sources are con- 
sistent with those expected for late-type galaxies, while the higher 
values for the radio galaxies suggest early-type hosts. A KS test 
finds a significant separation between the two samples, yielding a 
significance level of only p — 0.01 that the two subsamples are 
drawn from the same distribution. It should be noted that constrain- 
ing source morphology from concentration index typically requires 
higher-resolution data than that used in this study. However, the rel- 
atively clean separation seen here between these two sub-samples 
suggests a difference in morphology completely consistent with 
that implied by the Sersic index distributions produced by the 2- 
dimensional modelling. 

The ability of the two-dimensional modelling to recover accu- 
rate morphological information from simulated data, and the con- 
sistency of these results with multiple non-parametric tests of ri/2 
and basic morphology, all support the conclusion that the 2 ~ 2 ra- 
dio galaxies are already relaxed ellipticals, while the sub-millimetre 
galaxies are much better described as massive discs. 



6.3 Non-parametric tests of ri/2 

The Petrosian radius and ri/2 were determined from a curve- 
of-growth analysis for comparison with the modelling results. 
Surface-brightness profiles for each object were derived by using 
the IRAF package PHOT to measure aperture magnitudes through 



6.5 Image stacking 

As described in Section 3, the exposure times adopted in this study 
were selected to achieve a signal-to-noise of 5 at twice the antici- 
pated ~ 10 kpc scalelength of the radio galaxies, in order to allow 
accurate two-dimensional modelling of galaxy profiles. To explore 
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5 4 

Sersic index 

Figure 1. The distribution of Sersic index n fertile sub-millimetre galaxies 
(upper panel) and radio galaxies (lower panel) as derived from the two- 
dimensional modelling of the _R'-band images. 




0.4 0.6 

Axial ratio 

Figure 3. The distiibution of axial ratios for sub-millimetre galaxies (above) 
and radio galaxies (below) as derived from the two-dimensional modelling 
of the A'-band images. The secure (i.e. radio-identified) sub-millimetre 
galaxies are shown with dark grey shading, while the less secure optical- 
infrared identifications are indicated by light grey shading. 




5 10 

Effective radii / kpc 

Figure 2. The distribution of the half-light radius, ri/2> for the sub- 
millimetre galaxies (upper panel) and radio galaxies (lower panel), as de- 
lived from the two-dimensional modelling of the K-band images. 



the extent to which extended low-surface brightness emission could 
still have remained undetected, the image data for both galaxies 
sub-classes were stacked to create images with effective exposure 
times of 33.5 hours on UKIRT for the radio galaxies, and 21.7 hours 
on Gemini for the sub-millimetre galaxies (=80 hours on UKIRT). 
All the individual galaxy images were centroided, scaled to the av- 
erage luminosity, aligned to the same position angle, and mean- 
combined using the IRAF package IMCOMBINE. The resulting 
images were then modelled using PSFs created from a scaled stack 
of all individual PSFs. The stacked images, and the results of the 
two-dimensional model fits are presented in Fig.[5]and Table[5] 

The luminosity profiles of these two stacked images are shown 
in Fig. [6] It is clear that there is a significant amount of extended 
emission around the radio galaxies which can be traced out to a ra- 
dius of ~ 5 arcsec, consistent with a large relaxed elliptical. In 
contrast, the profile of the sub-millimetre galaxy stack confirms 
the less extended, power-law luminosity profiles of these objects, 
with no signal detected beyond a radius of 2 arcsec. This compar- 




Concentration index 

Figure 4. The distribution of concentration index for sub-millimetre galax- 
ies (above) and radio galaxies (below) as derived from the two-dimensional 
modelling of the A'-band images. The secure (i.e. radio-identified) sub- 
millimetre galaxies are shown with dark grey shading, while the less secure 
optical-infrared identifications are indicated by light grey shading. 



ison demonstrates very clearly that the NIRI imaging of the sub- 
millimetre galaxies is deep enough to reveal extended emission 
as displayed by the radio galaxies if it was present at a compara- 
ble level. The parameter values derived from the two-dimensional 
modelling of these two image stacks are given in Table |5] The re- 
sults are reassuringly consistent with the average values of the in- 
dividual object sub-samples which have been incorporated into the 
stacks, reinforcing the evidence that these two classes of galaxy are 
very different in terms of size and basic morphology. 

6.6 Analysis of degraded HST ACS /-band imaging 

The absence of significant morphological asymmetries in the K- 
band sub-millimetre galaxy modelling seems at first glance to be 
in conflict with the known presence of asymmetries seen in high- 



© 0000 RAS, IVINRAS 000, 000-000 



Properties ofsubmm and radio galaxies at high redshift 9 



t 



• 




Figure 5. The stacked ii'-band images of the radio galaxies (above) and sub-millimetre galaxies (below), along with the results of fitting these stacked images 
with the axi-symmetric two-dimensional galaxy models. The left-hand panels show the stacked images. The middle panels show the best-fitting models. The 
light-hand panels show the residual images produced by subtracting the best-fitting models from the stacked data. All panels are 8.0" X 8.0", and the images 
are shown with a linear greyscale in which black corresponds to 2. Sit above, and white to la below the median sky value in the stacked images. 
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Figure 6. The luminosity profiles (data and model fits) derived from the stacked ft'-band images of the z ~ 2 radio (left) and sub-millimetre galaxies (right). 
Model profiles are plotted out to the radius at which the signaknoise level in the stacked images allows them to be constrained. 
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Table 5. The results of the two-dimensional modelling of the stacked radio 
and sub-millimetre galaxy K-band images. 



Stack 



''1/2 < Me > Sersic n A"— mag 



At/ 



Radio galaxy 6.12 19.18 4.6 17.88 .97 

Sub-mm galaxy 2.19 18.65 1.4 19.48 .94 



resolution space-based HST imaging of sub-millimetre galaxies 
(e.g. Tacconi et al. 2008; Swinbank et al. 2010). The Swinbank 
et al. (2010) data reveals a population of often multi-component 
sources, typically not well described by standard "Hubble tuning 
fork" morphologies. One interpretation is, of course, that the galax- 
ies do indeed appear less disturbed when viewed at Xobs — 2.2 ^m 
than at shorter wavelengths. However, an obvious concern is that 
such strong asymmetries are in fact still present, but are simply not 
apparent in our 7f -band data due to the limitations of our (albeit 
good) ground-based seeing (typically ~ 0.5 arcsec FWHM). 

To attempt to distinguish between these alternatives, we se- 
lected the four objects in our ELAIS N2 sub-millimetre galaxy 
sample which were imaged at 7-band with HST ACS (PI: Almaini, 
PID: 9761), and which yielded optical detections of sufficient sig- 
nificance to warrant image modelling (namely ELAIS N2 850.1, 
N2 850.4, N2 850.7, and N2 850.8). We degraded these HST data 
to the same pixel-scale and seeing quality as our ground based ob- 
servations by convolving the ACS images with an averaged PSF 
from our ii'-band imaging programme. We then performed two- 
dimensional modelling of these "pseudo ground-based" red optical 
images in the same way as described above for the real A'-band 
imaging. 

These degraded images, and the results of trying to model 
them with axi-symmetric galaxy models are shown in Appendix 
A, Fig. A3. Unlike the A-band images, we found that the degraded 
ACS images of these four objects could not be so well-described by 
a simple axi-symmetric galaxy model, with much more significant 
asymmetries present in the model-subtracted residual images (with 
only 60 — 90% of the light accounted for by the axi-symmetric 
models). This favours the morphological k-correction interpreta- 
tion, supporting the view that the galaxies really are smoother when 
observed at A-band, longward of the rest-frame 4000A or B aimer 
break. This conclusion is consistent with the results of Swinbank 
et al (2010), who found that the /-band ACS images of their sub- 
millimetre galaxies generally displayed larger degrees of asymme- 
try than suggested by the NICMOS data. 

In conclusion, our modelling of the available data suggests 
that, while at least some z ~ 2 sub-millimetre galaxies may appear 
complex when viewed at rest-frame ultra-violet wavelengths (per- 
haps due to a mixture of patchy obscuration and/or localised star- 
burst activity), when viewed at 2.2 ^m, a single dominant galaxy is 
generally seen, the vast majority (> 95% in every case) of whose 
observed rest-frame optical emission can be well described by a 
simple, axi-symmetric exponential disc. 



7 DISCUSSION 

7.1 Galaxy morphologies 

One of the cleanest results from this study is that the radio galaxies, 
even at 2; ~ 2, appear to rather well described as de Vaucouleurs 
spheroids (i.e. Sersic index n = 4) whereas, in contrast, the sub- 
millimetre galaxies are much better described as exponential discs 



(i.e. Sersic index n = 1). This does not of course preclude the pos- 
sibility that the sub-millimetre galaxies could evolve into present- 
day ellipticals, but it is clear that, when selected as bright sub- 
millimetre emitters, the luminosity profiles of these galaxies are 
not well described by a de Vaucouleurs law (either individually, or 
in the stacked image shown in Fig. 5). The axial ratio distributions 
displayed by the two samples are consistent with this morphologi- 
cal division, although they are not formally statistically different. 

It is often argued that the extreme starbursts required to pro- 
duce easily detected sub-millimetre galaxies are the result of major 
mergers. For example, in the simulations of Narayan et al. (2009), 
the merging of two galaxies with stellar masses ~ 3 x 10^^ M0 
is required to produce a sub-millimetre flux density Ssso — 5 mjy. 
In addition, direct observational evidence that at least some sub- 
millimetre galaxies display complex, disturbed gas motions has 
been presented by Tacconi et al. (2008). 

However, it is not yet clear that major mergers are the domi- 
nant cause of high-redshift sub-millimetre galaxies. At least some 
sub-millimetre galaxies in fact display CO velocity gradients con- 
sistent with a compact, rotating disc (Tacconi et al. 2008), and 
Dave et al. (2010) have argued that the sub-millimetre galaxies 
predicted by current cosmological hydrodynamic simulations are 
massive galaxies being fed by smooth infall and gas-rich satellites 
at rates comparable to their star-formation rates. In this scenario the 
fraction of sub-millimetre galaxies expected to show clear evidence 
of a major merger is predicted to be small, and most sub-millimetre 
galaxies would be expected to display relatively undisturbed mor- 
phologies. 

From Fig. A2 it can be seen that while some of the sub- 
millimetre galaxies imaged here show evidence of a recent galaxy- 
galaxy interaction (e.g. N2850.4, N2850.8, LE850.3, LE850.7), in 
all cases the vast majority (>95%) of the A-band light from these 
galaxies is well described by a single, simple axi-symmetric galaxy 
model (see final column in Table 2 for individual object results). 

The average, deboosted 850 pan flux density of the sub- 
millimetre galaxies imaged here is (Sgso) — 6 miy. Our results in- 
dicate that, at least at this flux-density level, 2 ~ 2 sub-millimetre 
galaxies are massive, star-forming disc galaxies which generally 
display little obvious evidence of a recent massive galaxy-galaxy 
merger. The concept that sub-millimetre galaxies are simply the 
high-mass tail of the general star-forming galaxy population at 
z ~ 2 — 3 gains at least qualitative support from the comparison of 
the K-hand magnitudes measured here for sub-millimetre galaxies 
with the K-band magnitudes measured for Lyman-break galaxies 
at comparable redshifts by Shapley et al. (2001), shown in Fig.|7] 
How high-redshift sub-millimetre galaxies compare to present-day 
star-forming discs is discussed further below. 

7.2 Galaxy sizes 

The sub-millimetre galaxies are fairly compact, with stellar half- 
light radii ~ 3kpc. As already mentioned for individual sources 
above, this typical size is comparable to the gas sizes of sub- 
millimetre galaxies measured by Tacconi et al. (2008) using 
IRAM/PdB (ri/2 — 2kpc), a result which is at least consistent 
with the idea that both the stars and molecular gas are co-located in 
these star-forming discs. 

Until recently, such small sizes inferred for massive galaxies 
(Mstars ~ 2 X 10^^ Mq; see Section 7.3 below) at high-redshift 
might have been regarded as unexpected. However, recent work has 
shown that other populations of massive galaxies at 2 = 1.5 — 3 
are generally similarly compact, with ri/2 ~ 2 — 3 kpc (e.g. Daddi 
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et al. 2005; Tmjillo et al. 2006, 2007; Zirm et al. 2007; Cimatti 
et al. 2008; van Dokkum et al. 2008), with both star-forming and 
quiescent galaxies apparently a factor of 3 — 5 times smaller than 
present-day counterparts of comparable mass (Toft et al. 2009, but 
see also Manccini et al. 2010, and Saracco, Longhetti & Gargiulo 
2010). The controversial question of how such compact galaxies 
might evolve into present-day massive ellipticals is discussed fur- 
ther in the next sub-section, but observationally the key point is 
that the sizes of the sub-millimetre galaxies measured here are not 
unexpectedly small in the context of those measured for other com- 
parably massive galaxies at comparable redshifts. 

The radio galaxies are markedly larger, with stellar half-light 
radii ~ 8kpc. However, they are also typically ~ 3 times more 
massive. Interestingly, the average size of the z > 1.5 radio galax- 
ies measured here is very similar to that displayed by the bright- 
est cluster galaxies at 2 ~ 1.5 recently studied by by Collins et 
al. (2009). Further morphological studies of purely mass-selected 
high-mass (Alstar > 5 x lO'^^ M©) galaxies at z ~ 2 are required 
to establish whether these larger scale sizes are typical of this ex- 
treme mass regime at these redshifts, or whether they are seen in 
only certain subclasses of object which may have special merger 
histories (radio selection, or brightest cluster galaxy selection may 
preferentially select the products of early major mergers). 

However, while clearly larger than the sub-millimetre galax- 
ies, the 2; ~ 2 radio galaxies are nonetheless undoubtedly more 
compact than comparably radio-luminous galaxies at 2 < 1. 
Specifically, the 2 ~ 2 radio galaxies studied here are found to 
be, on average, a factor ~ 1.5 smaller than the < 2 < 1 3CR 
radio galaxies studied by McLure et al (2004). Motivated by the 
fact that the 2 ~ 2 radio galaxies appear to be well described as 
dynamically relaxed spheroids, we have attempted to explore fur- 
ther whether they follow the scaling relations displayed by ellip- 
tical galaxies in general, and radio galaxies in particular, at lower 
redshifts. 



7.2.1 The Kormendy relation 

Early-type galaxies are known to exist on a two-dimensional man- 
ifold (the "fundamental plane"), defined by effective scalelength, 
mean surface brightness, and the central stellar velocity dispersion 
in three-dimensional parameter space (e.g. Dressier et al. 1987; 
Djorgovski & Davis 1987). In the absence of velocity dispersion 
measurements, the Kormendy or /ie — r^ relation (Kormendy 1977) 
allows examination of the photometric projection of the fundamen- 
tal plane using only the scalelength and surface-brightness parame- 
ters. The i?— band Kormendy relation deduced for the 2 ~ 2 radio 
galaxies studied here is compared with that deduced for 2 = 
3CR-type radio galaxies by McLure et al. (2004) in the upper panel 
of Fig. [8] The < /i >e values of the 2 ~ 2 radio galaxies have 
been corrected for cosmological surface-brightness dimming, and 
fc-corrected into the _R— band. The data points from the McLure et 
al. sample have been corrected to 2 ~ assuming passive evo- 
lution derived from the Bruzual & Chariot (2003) galaxy-evolution 
models. The best-fitting form of the Kormendy relation for the non- 
evolutionary corrected 2 ~ 2 radio galaxy sample is 



< ^l >e= 3.96(±1.19) log re + 14.39(±1.03) 



(2) 



To within the errors, the slope of this relation is indistinguish- 
able from that displayed by the 2 ~ radio galaxy sample. If the 
two samples are fitted with a Kormendy relation of fixed slope 3.87 
(intermediate between the slopes of their two independent fits), 
as shown in the upper panel of Fig. [8] then the relations appear 



to be offset by a (vertical) shift in rest-frame _R-band luminos- 
ity equivalent to 1.76 magnitudes. This is in very good agreement 
with the expected passive evolutionary correction as derived from 
a zform ~ 3 burst model, using the evolutionary synthesis models 
of Bruzual & Chariot (2003) (A_R ~ 1.65). Motivated by this, we 
have then applied a passive evolution correction, galaxy-by-galaxy, 
to the 2 ~ 2 radio-galaxy sample, to move them to 2 ~ assum- 
ing no evolution in scalelength. The results are over-plotted on the 
2 ~ data in the lower panel of Fig.jS] and a fit to the combined 
dataset yields 



< fi >e= 3.74(±0.44) log re + 16.38(±0.47) 



(3) 



In summary, not only are the 2 ~ 2 radio galaxies well- 
described by a de Vaucouleurs r^'* law, but they also display a 
Kormendy relation which is indistinguishable from that displayed 
by low-redshift radio galaxies (and other radio-quiet massive ellip- 
ticals) after simply correcting their luminosities for anticipated pas- 
sive evolution of their stellar populations. This suggests that any 
further increase in the typical sizes of these objects since 2 ~ 2, 
must be accompanied by sufficient increase in mass to preserve the 
form of the Kormendy relation. 

In concluding this discussion of the radio galaxy sizes and 
morphologies, it is interesting to assess how our results compare 
with those of Pentericci et al. (2001), who attempted to fit models 
to 14 of the 19 radio galaxies in the redshift range 2 — 1.7 — 3.0 
which they imaged with HST NICMOS in the //-band. They re- 
ported that this did not, in general, work well and that many of the 
highest-redshift radio galaxies had complex, clumpy morphologies. 
But they did successfully fit de Vaucouleurs models to 5 of their 
objects, all of which lay at the low-redshift end of their sample, 
at 2 ~ 1.8 — 2.2 (more comparable to the redshifts of the radio 
galaxies studied here). The average size of these five radio galaxies 
was found to be 5.5 kpc, again consistent with our findings in the 
present study. 

These results, now reinforced by our own findings, lend sup- 
port to the general idea that redshift 2 ~ 2 corresponds to the 
epoch when the most massive galaxies have evolved into relaxed 
spheroids. However, we note that it would be interesting to extend 
the ultra-deep A'-band imaging of radio galaxies out to 2 ~ 3, to 
explore the extent to which the NICMOS imaging of the higher red- 
shift end of the Pentericci et al. sample could have been affected by 
emission-line contamination, the movement of the 4000A/Balmer 
break through the F160W filter, and the rapidly decreasing sensi- 
tivity to low-surface brightness emission. 



7.3 Surface mass densities 

Finally we use our A'-band galaxy model fits to explore where the 
2 ~ 2 radio galaxies and sub-millimetre galaxies lie in the stellar 
surface-mass-density versus stellar mass plane. 

The stellar masses (Alstar) of the galaxies were estimated by 
matching the measured total (i.e. 6-arcsec diameter aperture, ex- 
cluding companions) if-band magnitudes of the galaxies to the 
predictions of the evolutionary synthesis models of Bruzual & 
Chariot (2003). Given the evidence that the radio galaxies are al- 
ready evolving essentially passively, a uniform burst-model age of 
2 Gyr was adopted in establishing their A'-band mass-to-light ra- 
tio. The inferred average stellar mass of the 2 ~ 2 galaxies is 
~ 7 X 10^^ M0, consistent with the findings of McLure et al. 
(2004). 

For the sub-millimetre galaxies, which are by definition not 
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Figure 7. The ii'-band magnitude distribution of Lyman-break galaxies 
at 2 ~ 2 — 3 (upper panel) from Shapley et al. (2001), compared with 
that derived here for sub-millimetre galaxies (lower panel) at comparable 
redshifts. The radio-identified sub-milUmetre galaxies are shown with dark 
grey shading, while the optical-infrared identifications are indicated by light 
grey shading. 

passively evolving, such a simple assumption is clearly not justi- 
fied. However, for the majority of these galaxies, multi-wavelength 
(optical-mid-infrared) photometry is now available from Subaru 
and Spitzer which, when combined with the i('-band photometry 
presented here, allows detailed fitting to the predictions of the evo- 
lutionary synthesis models (e.g. Dye et al. 2008). The stellar masses 
derived from these fits were used to calibrate the typical mass- 
to-light ratio for the sub-millimetre galaxy sample. This was then 
used to estimate the masses of the remaining sub-millimetre galax- 
ies for which less complete multi-frequency photometry was avail- 
able. The inferred average stellar mass of the z ~ 2 sub-millimetre 
galaxies is ~ 2.2 x 10^^ Mq, consistent with the results of Dye et 
al. (2008) and Schael et al. (2010). 

The average stellar mass surface density, ctsq, is then defined 
as 

0.5Mstar ,., 



O"50 ~ 



'1/2 



where ri/2 is the half-light radius of the starlight, as derived 
from our axi-symmetric two-dimensional fits to the iC-band images 
(see Table 2). 

The results are plotted in Fig.|9] following Zirm et al. (2007), 
alongside comparable measurements made for other high-redshift 
populations taken from a number of recent studies (Daddi et al. 
2005, Longhetti et al. 2007, Zirm et al. 2007, van Dokkum et al. 
2008). 

As discussed by Zirm et al. (2007), the star-forming DRGs 
and Lyman-break galaxies at z ~ 2 appear to lie within within 
the locus defined by local galaxies, as indicated by the grey shaded 
region in Fig. |9](Shen et al. 2003). However, the more quiescent 
high-redshift galaxies appear to lie well above this relation. This 
is simply a graphical restatement of the fact that, given their in- 
ferred stellar masses, these objects are surprisingly compact com- 
pared to low-redshift galaxies. The question of how such objects 
can have evolved onto the local galaxy locus has been considered 
by several authors. For example, it has been suggested that while 
secular evolution would exceed the Hubble time, the necessary re- 
laxation of ri/2 could arise via dissipationless 'dry' merging (e.g. 
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Figure 8. The top panel shows the best-fitting fixed-slope Kormendy rela- 
tions for the 2 ~ 2 radio galaxy sample (solid line) and the z ~ 3CR-type 
galaxies from McLure et al. (2004) (dashed line). In the bottom panel the 
2 ~ 2 radio galaxies have been corrected for passive evolution of their stel- 
lar populations between z ~ 2 and z ~ as described in the text, and 
the solid line indicates the best-fitting Kormendy relation to the combined 
sample. 



van Dokkum 2005), with a linear decrease in the surface density 
with accumulated mass (Nipoti, Londrillo, & Ciotti 2003). A fidu- 
cial normalisation of this power-law trend is shown in Fig. |9] the 
implication is that dynamical evolution via dry mergers would al- 
low the quiescent DRGs to reach local galaxy densities, but only at 
very high masses (~ 10^'^ Mq). 

What can be learned about the z ~ 2 radio galaxies and sub- 
millimetre galaxies from their inferred location on this diagram? 
First, we note that most (but not all) of the z ~ 2 radio galax- 
ies already lie on the local galaxy locus, with project stellar mass 
densities of ctso — 10® MqIcpc"'^. This is interesting, but perhaps 
not unexpected given their moderately large sizes and apparently 
dynamically relaxed stellar populations. Four of the radio galaxies 
have higher surface mass densities. This might be telling us that 
their stellar masses have been over-estimated, or their half-light 
radii under-estimated, but the two obvious cases of an AGN contri- 
bution have already been removed and are not plotted in Fig.|9] Al- 
ternatively these objects could really be destined to evolve further, 
in which case the dry-merging tracks shown in Fig. |9] would im- 
ply that they will land on the present day locus with stellar masses 
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Figure 9. Average surface mass density within the half-light radius plotted versus stellai' mass. The z ~ 2 radio galaxies studied here are indicated by the 
filled circles, while the sub-millimetre galaxies are indicated by the solid stars (black for radio-identified galaxies, grey for optical-infrared identifications). 
Open circles indicate the location of passively evolving massive early-type galaxies in the redshift range 1.4 < z < 2.6 from Daddi et al. (2005), Longhetti 
et al. (2007), Zirm et al. (2007), and van Dokkum et al. (2008) are shown. The measurements for star-forming distant red galaxies (DRGs) and Lyman-break 
galaxies given by Zirm et al. (2007) are also plotted as large open stars and small open stars respectively. The solid line traces the projected trend of surface 
density with mass expected for dissipationless mergers (Nipoti, Londrillo, & Ciotti 2003) normalised to the quiescent DRGs from Zirm et al. (2007). The 
dotted lines show the same trend if the stellar mass estimates from Zirm et al. (2007) are systematically too high or low by a factor of 3. Light-grey shading 
represents local values for early- and late-type galaxies from Shen et al. (2003). 



~ 2 — 3 X 10 Mq . This is in fact not unreasonable, given that 
they are already among the most massive galaxies known at z ~ 2, 
and yet still appear somewhat more compact than the most massive 
elliptical galaxies in the present-day universe. 

Interestingly, and unlike the lower-luminosity star-forming 
Lyman-break galaxies, most of the sub-millimetre galaxies also 
appear to lie above the local galaxy locus, in the same general 
regime as the moderately compact quiescent galaxies. Thus, while 
the morphological evidence suggests the sub-millimetre galaxies 
are star-forming discs, in terms of surface mass density they can 
clearly evolve into high-redshift quiescent galaxies without sig- 
nificant size evolution. Also, if they subsequently follow the dry- 
merger evolutionary tracks indicated in Fig. |9l they are ultimately 
destined to evolve into present-day massive elliptical galaxies with 

Mstar > 5 X 10" Mq. 

There are 5 sub-millimetre galaxies which appear to have 
lower surface mass densities, and thus lie within the locus defined 
by present-day galaxies and high-redshift DRGs and Lyman-break 
galaxies. Both of the radio unidentified sub-millimetre galaxies lie 
within this subsample, suggesting either that they have been mis- 
identified, or (not unreasonably) that strong radio emission is an 
indicator for the densest starbursts. 

In summary. Fig. |9] supports the conclusion that, while the 
sub-millimetre galaxies appear (at the epoch of observation) to be 
best described as star-forming discs, their stellar masses and stel- 



lar mass surface densities indicate that they are capable of evolving 
rapidly into dense, quiescent galaxies at z > 1.5, and ultimately 
into present-day massive ellipticals. 



8 CONCLUSIONS 

We have obtained deep, high-quality A'-band images of complete 
subsamples of powerful radio and sub-millimetre galaxies at red- 
shifts 2 ~ 2. The data were obtained in the very best available see- 
ing (FWHM ~ 0.5 arcsec) through the queue-based observing sys- 
tems at the United Kingdom InfraRed Telescope and Gemini North, 
with integration times scaled to ensure that comparable rest-frame 
surface brightness levels are reached for all of the galaxies. We 
have analysed the resulting images by fitting two-dimensional axi- 
symmetric galaxy models to determine basic galaxy morphologi- 
cal parameters at rest-frame optical wavelengths \rest > 4000A, 
varying luminosity, half-light radius ri/2, and Sersic index n. We 
have also undertaken an extensive set of simulations which show 
that, for galaxies with scalelengths ri/2 > 2 kpc, the modelling of 
our deep, high-quality ground-based images should yield half-light 
radii and Sersic indices accurate to ~ 10%. 

We find that, while a minority of the images show some evi- 
dence of galaxy interactions, > 95% of the rest-frame optical light 
in all the galaxies is well-described by these simple axi-symmetric 
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models. We also find evidence for a clear difference in morphology 
between the radio galaxies and the sub-millimetre galaxies. Specif- 
ically, the fits to the individual images reveal that the radio galaxies 
are moderately large ((ri/2) = 8.4± 1.1 kpc; median r-i/2 ~ 7.8), 
de Vaucouleurs spheroids ((n) — 4.07 ± 0.27; median n = 3.87), 
while the sub-millimetre galaxies appear to be moderately compact 
(('"1/2) ~ 3.4 ± 0.3 kpc; median r-i/2 ~ 3.1 kpc) exponential discs 
((n) = 1.44 ± 0.16; median n = 1.08). Model fits to the stacked 
images yield very similar results. 

We find that the z ~ 2 radio galaxies display a well-defined 
Kormendy (fie — re) relation, which is offset in surface brightness 
from the local relation by ~ 1.7 mag. This luminosity offset is 
consistent with that expected due to purely passive evolution be- 
tween z ~ 2 and z ~ for a stellar population formed at z > 3. 
However, while larger than other recently-studied massive galaxy 
populations at comparable redshifts, the z ~ 2 radio galaxies are 
still on average a factor ~ 2 times smaller than their local counter- 
parts, suggesting that they have still to undergo further dynamical 
evolution. 

The scalelengths we have derived for the starlight in the 
sub-millimetre galaxies are consistent with those reported for the 
molecular gas (e.g. Tacconi et al. 2008), suggesting that the two 
may be co-located. Such scalelengths are also similar to those 
of comparably-massive quiescent galaxies at z > 1.5, allow- 
ing the possibility of an evolutionary connection following cessa- 
tion/quenching of the observed star-formation activity. 

Finally, in terms of stellar mass surface density, we find that 
the majority of the radio galaxies lie within the locus defined by 
local ellipticals of comparable stellar mass. In contrast, while best 
modelled as discs at the epoch of bright dust-enshrouded star for- 
mation, most of the sub-millimetre galaxies have higher stellar 
mass densities than local galaxies, and appear destined to evolve 
into present-day massive ellipticals. 

With the advent of Herschel and SCUBA-2, coupled with the 
wide-field near-infrared imaging provided by VISTA and WFC3/IR 
on HST, it should soon be possible to extend this type of study to 
larger samples of sub-millimetre galaxies which include even more 
extreme starbursts at z ~ 2. It will then be possible to test whether 
the dense disc morphologies of sub-millimetre galaxies uncovered 
in the present study are a universal feature of dust-enshrouded star- 
bursts at high redshift, or whether they are specific to the mass/star- 
formation rate regime sampled in the early sub-millimetre surveys. 
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APPENDIX A: SCIENCE IMAGES, GALAXY MODEL 
FITS, AND RESIDUALS 

In this appendix we present greyscale plots of the ii'-band im- 
ages, the axi-symmetric model fits, and then the residual images 
after subtraction of the PSF-convolved model. Plots for the radio 
galaxies are shown in Fig. Al, with the plots for the sub-millimetre 
galaxies presented in Fig. A2. 

In addition, in Fig. A3, we show comparable greyscale plots 
of the HST ACS /-band images of the sub-milimeater galaxies 
N2850.1, N2850.4, N2850.7, and N2850.8, again accompanied by 
the axi-symmetric model fits, and the model-subtracted residual im- 
ages. 
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Figure Al. Two-dimensional modelling of the z ~ 2 radio galaxies. The left-hand panel shows the ii'-band image centred on the radio galaxy. The middle 
panel shows the best-fitting two-dimensional model. The right-hand panel shows the residual image after subtraction of the model from the data. All panels 
are 12.0" x 12.0", and the images are shown with a linear greyscale in which black corresponds to 2.5a- above, and white to Icr below the median sky value. 
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Figure Al. - continued 
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Figure Al. - continued 
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Figure A2. Two-dimensional modelling of the sub-millimetre galaxies. The left-hand panel shows the ii'-band image centred on the sub-millimetre galaxy. 
The middle panel shows the best-fitting two-dimensional model. The right-hand panel shows the residual image after subtraction of the model from the data. 
All panels are 12.0" X 12.0", and the images are shown with a linear greyscale in which black corresponds to 2.5a" above, and white to la below the median 
sky value. 
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Figure A2. - continued 
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Figure A2. - continued 
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Figure A2. - continued 
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Figure A3. Two-dimensional modelling of tlie HST ACS 7-band images of ELAIS N2 850. 1 , 850.4, 850.7 and 850.8 (PI Almaini, PID: 9761), after degradation 
to ground-based seeing of ~ 0.5 arcsec (FWHM). The left-hand panel shows the smoothed ACS image re-sampled to Gemini pixel-scale. The middle panel 
shows the best-fitting two-dimensional model. The right-hand panel shows the residual image after subtraction of the model from the data. All panels are 12.0" 
X 12.0", and the images are shown with a linear greyscale in which black corresponds to 2. Sir above, and white to Ict below the median sky value. 
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